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Scaling  Laws  for  Diamond  Chemical  Vapor  Deposition. 
II:  Atomic  Hydrogen  Transport 


D.  G.  Goodwin 

Division  of  Engineering  and  Applied  Science 
California  Institute  of  Technology 
Pasadena,  CA  91125 


Scaling  relations  are  developed  to  allow  estimation  of  the  atomic  hydro¬ 
gen  concentration  at  the  substrate  during  diamond  chemical  vapor  deposition  for 
both  diffusion-dominated  and  convection-dominated  reactors.  In  the  convection- 
dominated  case,  it  is  shown  that  there  exists  an  optimal  Mach  number  which  maxi¬ 
mizes  the  H  concentration  delivered  to  the  substrate.  In  addition,  when  homogeneous 
recombination  is  taken  into  account,  there  exists  an  optimal  operating  pressure.  This 
analysis  shows  that  a  sonic  flow  of  highly  dissociated  hydrogen  at  a  pressure  near 
atmospheric  is  optimal  for  rapid  growth  of  high-quality  diamond. 


I.  INTRODUCTION 

To  scale-up  diamond  chemical  vapor  deposition  processes,  two  factors  must  be 
considered.  The  dependence  of  film  growth  rate  and  quality  on  the  local  chemical 
environment  at  the  substrate  must  be  known,  and  the  relation  of  this  environment  to 
controllable  process  parameters  must  also  be  understood.  In  a  companion  paper  to 
this  one  [1],  we  have  presented  a  simplified  model  of  diamond  surface  chemistry,  which 
results  in  simple  approximate  expressions  for  the  growth  rate  and  defect  density  in 
terms  of  the  concentrations  of  H  and  CH3  at  the  substrate.  It  was  shown  that  a  simple 
defect  formation  model,  in  which  defects  are  assumed  to  result  from  interaction  of 
surface  species,  leads  to  a  relationship  between  the  defect  density  Xdej,  the  linear 
growth  rate  G,  and  the  atomic  hydrogen  concentration  at  the  surface  (H]o: 


where  n  «  2.  This  expression  shows  that  increasing  the  growth  rate  (for  example, 
by  adding  additional  hydrocarbon)  requires  a  compensating  increase  in  the  atomic 
hydrogen  concentration  to  maintain  the  same  film  quality.  This  qualitative  conclusion 
is  in  accord  with  much  experience  in  diamond  film  growth,  and  with  a  similar  analysis 
by  Butler  and  Woodin  [2]. 

Therefore,  a  critical  issue  for  process  scale-up  and  optimization  is  to  determine 
the  scaling  of  [Hjo  with  process  parameters.  This  question  is  complex,  since  the  H 
concentration  at  the  substrate  is  determined  by  a  balance  between  surface  recombi¬ 
nation  and  transport  to  the  surface,  either  by  diffusion  or  convection,  and  may  be 
affected  by  homogeneous  recombination  as  well. 

In  this  paper,  the  question  of  H  transport  under  diamond  chemical  vapor  depo¬ 
sition  (CVD)  conditions  is  examined  in  detail.  Simple  scaling  relations  are  derived 
for  the  effects  of  pressure,  flow  velocity,  substrate  diameter,  and  degree  of  gas  ac¬ 
tivation  on  the  H  concentration  at  the  substrate,  for  both  diffusion-dominated  and 
convection-dominated  reactors. 

After  first  deriving  the  basic  relations  valid  for  all  reactors,  transport  in  low- 
pressure  diffusion-dominated  reactors  is  considered,  followed  by  a  detailed  discussion 
of  H  transport  in  high-growth-rate,  convection-dominated  reactors  (primarily  plasma 
and  combustion  torches).  In  the  latter  case,  it  is  shown  that  there  exists  an  optimal 
Mach  number  and  pressure  for  H  transport.  Expressions  are  also  given  for  the  heat 
flux  to  the  substrate. 

For  most  of  the  results  presented  here,  simple  transport  theory  is  used  to  derive 
the  scaling  relations.  In  some  cases,  particularly  when  gas-phase  chemistry  is  con¬ 
sidered,  these  results  are  supplemented  by  numerical  simulations  of  the  flew.  In  this 
way,  we  may  assess  the  importance  of  effects  (such  as  variable  transport  properties, 
large  temperature  gradients,  and  gas-phase  chemistry)  which  we  neglect  in  the  simple 
theoretical  derivations. 

The  code  used  for  the  simulations  is  one  we  have  developed  and  used  extensively 
for  diamond  CVD  modeling  [3-7].  In  this  code,  the  governing  transport  equations 
for  axisymmetric  flow  are  solved  with  appropriate  boundary  conditions  to  determine 
the  velocity,  temperature,  and  concentration  profiles  near  the  surface.  Both  het¬ 
erogeneous  and  homogeneous  chemistry  are  included,  along  with  variable  transport 
properties  calculated  using  the  subroutines  of  Kee  et  al.  [8].  Green  et  al.  [7]  have 
recently  compared  species  and  temperature  profiles  in  the  stagnation-point  boundary 
layer  of  an  atmospheric-pressure,  inductively-coupled  plasma  torch  calculated  using 
this  code  with  measured  profiles  obtained  with  degenerate  four-wave  mixing,  and  find 
good  agreement. 

The  general  equations  derived  here  are  evaluated  numerically  for  the  case  of  a 
dilute  hydrocarbon  in  hydrogen  gas  mixture.  These  evaluated  formulas  are  given  in 
mixed  units,  with  lengths  in  cm  and  the  pressure  in  Torr.  The  substrate  temperature 
is  assumed  to  be  1200  K  in  all  cases. 
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II,  BASIC  EQUATIONS 


For  the  purposes  of  species  transport  to  the  surface,  the  Peclet  number  provides  a 
useful  means  of  classifying  different  diamond  growth  methods  [9],  The  Peclet  number 
is  defined  by  Pe  =  uL/D,  where  u  is  a  characteristic  flow  velocity,  L  a  characteristic 
length  scale  (for  example,  the  substrate  diameter),  and  D  is  the  diffusion  coefficient 
for  the  species  in  question.  If  Pe  <|C  1,  species  transport  is  dominated  by  diffusion, 
and  is  independent  of  the  gas  velocity.  Low-pressure,  low-velocity  diamond  CVD 
reactors  (for  example,  conventional  hot  filament  and  microwave  systems)  operate  in 
this  regime  (9, 3, 10, 11]. 

In  the  opposite  limit  Pe  >•  1,  convective  species  transport  is  dominant.  High- 
growth-rate  methods  (combustion  torches,  plasma  jets)  typically  operate  in  this 
regime.  In  this  case,  transport  of  species  in  regions  far  from  any  surface  is  primarily 
due  to  bulk  fluid  motion.  However,  even  if  Pe  1  there  is  always  a  thin  fluid  layer 
near  all  solid  surfaces  where  diffusive  transport  dominates,  since  the  velocity  must  go 
to  zero  at  the  surface.  This  thin  layer,  through  which  species  must  diffuse  to  reach  the 
surface,  is  the  boundary  layer.  The  characteristics  of  the  boundary  layer  determine 
H  transport  to  the  substrate  in  these  systems. 

A  parameter  similar  to  the  Peclet  number  is  the  Reynolds  number,  defined  by 
Re  =  uL/i/,  where  v  is  the  kinematic  viscosity  of  the  gas.  The  Reynolds  number 
determines  the  nature  of  the  velocity  field,  and  governs  the  transition  to  turbulence. 
For  gases,  Re  and  Pe  are  similar  in  magnitude,  since  the  ratio  i //£>  (the  Schmidt 
number)  is  of  order  one.  For  example,  for  H2  diffusing  through  air,  ujD  —  0.22,  while 
for  C02  diffusing  through  H2  vj D  =  1.58  [12].  A  third  transport  property  important 
for  heat  transfer  is  the  thermal  diffusivity  a.  This,  too,  is  similar  in  magnitude  to  v 
for  most  gases.  For  most  diatomic  gases  at  reasonable  temperatures,  the  ratio  vja 
(the  Prandtl  number)  is  approximately  0.7. 

Since  the  gas  velocity  is  zero  at  the  surface  whatever  the  value  of  Pe  (neglecting 
the  small  Stefan  velocity  due  to  net  mass  deposition),  the  molar  flux  of  H  to  the 
surface  jjj  may  be  expressed  as  [13] 

JH  =  n°DH{dir)M  (2) 

where  n0  =  p/RTq  is  the  molar  gas  density  at  the  surface,  Dh  is  the  diffusion  coef¬ 
ficient  for  atomic  hydrogen,  Xu  is  the  atomic  hydrogen  mole  fraction,  and  2  is  the 
distance  from  the  surface. 

At  steady-state,  the  flux  of  H  to  the  surface  will  just  balance  the  loss  of  H  due  to 
surface  recombination,  and  therefore 

3h  =  Rhi  (3) 

where  Rh  is  the  surface  recombination  rate,  given  by 
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(4) 


Rh  =  7  H 


[H]cw 
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Here  7 h  is  the  recombination  coefficient,  and  c«  is  the  mean  thermal  speed  of  an  H 
atom. 

Recently,  two  groups  have  reported  measurements  of  7//  at  temperatures  near 
1200  K  [14, 15].  Both  studies  indicate  7 h  at  this  temperature  is  slightly  greater  than 
0.1,  indicating  that  H  recombines  readily  on  the  diamond  surface. 

Substituting  from  Equations  (2)  and  (4)  into  Eq.  (3),  the  H  mole  fraction  at  the 
surface  Xh,o  may  be  expressed  as 


Xh,  0 


4  Dh  f  dXu  \ 

1HCH  \  dz  )  z=Q  ' 


The  H  concentration  at  the  surface  is  then  given  by 


[H]0  =  tioXh,o- 


(5) 

(6) 


Let  us  define  a  characteristic  diffusion  length  scale  id  such  that 


XH,rej  ~  Xfjfi 


(7) 


where  XnyTtj  is  the  H  mole  fraction  at  some  specified  “reference”  position  (for  exam¬ 
ple,  the  filament  surface  in  a  hot  filament  reactor,  or  a  point  outside  the  boundary 
layer  in  a  high-Pe  reactor).  Then  Eq.  (2)  may  be  rewritten  as 


Jh  = 


«o  Dh 

U 


{Xh, ref  ~  Xni0), 


(8) 


and  Eq.  (5)  as 

(9) 

The  quantity  £j  is  clearly  an  important  parameter  for  H  transport:  H  at  the 
surface  will  be  maximized  when  id  is  made  as  small  as  possible.  (In  the  mass  transfer 
literature,  the  nondimensional  ratio  where  L  is  a  characteristic  length  scale  for 
the  problem,  is  known  as  the  “Sherwood  number”  [12, 13].)  Note  that  Eq.  (9)  is  always 
valid,  since  it  is  based  only  on  the  definition  of  id-  The  value  of  id ,  or  equivalently,  the 
Sherwood  number,  depends  on  all  properties  of  the  flow  (geometry,  Peclet  number, 
etc.)  and  is  determined  by  solving  the  governing  transport  equations.  If  Pe  <C  1  and 
homogeneous  recombination  can  be  neglected,  then  id  depends  primarily  on  reactor 
geometry;  on  the  other  hand,  if  Pe  >  1 ,  then  the  flow  velocity  and  gas  properties  also 
affect  id- 


Xh, 0  _  ^  7 HidCH 

Xh,™/  $Dh  . 
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Note  that  from  Eq.  (8),  the  quantity  tj/noDu  may  be  interpreted  as  a  “resistance” 
for  H  transport,  with  the  mole  fraction  difference  as  the  driving  potential.  Since  n0Dn 
is  independent  of  pressure  and  only  weakly  dependent  on  temperature,  the  resistance 
for  H  transport  is  a  function  primarily  of  Ia.  Explicit  values  for  £a  will  be  given  below. 

Let  us  define  another  length  scale  A//  as 


ch 


(10) 


This  length  scale  is  a  function  only  of  temperature,  pressure,  and  gas  composition, 
and  is  independent  of  flow  velocity.  Since  from  simple  kinetic  theory  Dh  is  roughly 
given  by  (l/3)cf/Am/p  [13],  where  Xmjp  is  the  mean  free  path,  A//  is  approximately 
equal  to  the  mean  free  path.  In  fact,  we  may  take  Eq.  (10)  to  define  the  mean  free 
path  for  H  transport.  Then  Eq.  (9)  may  be  written  as 

Xh,o 


X H,ref 


1  + 


7  H 
KnJ 


-l 


where 


K..1S 


(11) 


(12) 


is  the  Knudsen  number  for  H  based  on  l a. 

In  the  limit  Kn  >•  7//,  Eq.  (li)  reduces  to  Xu, a  =  Xh,tCj-  In  this  limit,  spatial 
gradients  in  H  mole  fraction  are  negligible.  In  the  opposite  limit  (Kn  7h),  the  H 
mole  fraction  at  the  surface  is  much  less  than  that  at  the  reference  location.  In  this 
limit,  gradients  in  H  mole  fraction  are  large,  and  Eq.  (11)  reduces  to 


Xh, 0  _  Kn 
-Xtf.re/  7  H 


(13) 


This  limit  defines  the  regime  in  which  H  transport  is  diffusion-limited. 

The  value  of  Xu  may  be  estimated  for  typical  CVD  diamond  growth  conditions. 
For  atomic  hydrogen  diffusing  through  a  bath  gas  of  H2,  Dh  (cm2/s)  at  1200  K 
is  approximately  1.7  x  10 4/p  [8],  where  p  is  the  pressure  in  Torr.  (A  few  percent 
hydrocarbon  mole  fraction  does  not  significantly  change  Dh-)  At  this  temperature, 
ch  =  5  x  105  cm/s;  therefore, 

0.14 

A h  ~ -  cm,  (14) 

V 

with  p  in  Torr.  Since  H  transport  begins  to  become  diffusion-limited  when  Kn  ss  7#, 
if  7 h  «  0.1,  then  at  20  Torr  Ia  must  be  less  than  0.7  mm  to  avoid  diffusion  limitations 
to  H  transport;  at  1  atm,  Ia  would  have  to  be  less  than  18  pm.  As  discussed  below, 
for  most  current  diamond  CVD  methods  Ia  is  not  this  small;  therefore,  in  most  cases 
H  transport  is  diffusion  limited,  and  we  would  expect  large  spatial  gradients  in  H 
concentration  within  the  reactor.  In  this  case,  to  increase  H  at  the  surface  £a  should 
be  made  smaller;  we  consider  below  means  to  do  this  both  for  low-Pe  and  high-Pe 
reactors. 
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III.  DIFFUSION-DOMINATED  METHODS 

If  Pe  <C  1,  then  id  is  determined  by  the  solution  of  the  diffusion  equation.  If  the 
diffusion  coefficient  is  constant  and  homogeneous  recombination  can  be  neglected, 
then  id  is  a  function  only  of  reactor  geometry.  Some  values  for  simple  geometries  are 
given  in  Table  I.  Note  that  for  cylinderical  and  spherical  geometries  Id  at  the  outer 
surface  can  be  significantly  greater  than  the  radius,  due  to  the  increase  of  area  with 
radius. 

It  will  be  shown  in  the  next  section  that  convection-dominated,  high-Pe  methods 
can  easily  achieve  4  C  1  mm.  Comparable  surface  H  concentrations  (and  therefore 
growth  rates)  could  in  principle  be  obtained  in  low-Pe  systems  if  the  geometry  were 
such  that  a  similar  id  resulted.  We  show  below  that  in  conventional  hot-filament  and 
microwave  plasma  systems,  id  is  of  order  several  centimeters,  although  by  placing  the 
source  of  atomic  hydrogen  closer  to  the  substrate,  smaller  values  could  be  obtained. 

A.  Hot-Filament  Reactors 

Let  us  estimate  id  and  the  H  mole  fraction  at  the  substrate  for  a  “typical”  hot- 
filament  reactor.  Quantitative  measurements  of  H  mole  fraction  in  this  environment 
have  been  made  by  Schafer  et  al.  [16]  and  Chen  et  al.  [17]  using  non-intrusive  optical 
methods,  and  by  Hsu  [18-20]  using  molecular-beam  mass  spectroscopy,  sampling  gas 
through  a  hole  in  the  substrate.  The  optical  measurements  show  that  the  H  concen¬ 
tration  decreases  with  increasing  distance  from  the  filament,  clearly  demonstrating 
diffusion-limited  H  transport  away  from  the  filament  [16, 17]. 

In  the  experiments  of  Hsu,  the  conditions  employed  were  a  pressure  of  20  Torr,  an 
uncoiled  filament  consisting  of  a  tungsten  wire  of  diameter  0.25  mm  and  length  5  cm, 
and  a  filament-substrate  separation  of  1.3  cm.  As  a  first  approximation,  let  us  assume 
that  the  filament  is  infinitely  long,  and  that  the  substrate  is  an  infinitely-long  cylinder 
of  radius  1.3  cm,  concentric  with  the  filament.  From  Table  I,  taking  R\  =  0.0125  cm 
and  /?2  =  1.3  cm,  we  find  £d  «  6  cm.  At  20  Torr,  we  then  find  from  Eq.  (10) 
that  A h  =  70  pm,  and  from  Eq.  (11)  we  estimate  Xn,o /Xhju  ~  0.011,  assuming 
7h  =  0.1.  If  it  is  assumed  that  the  H  mole  fraction  at  the  filament  corresponds  to 
chemical  equilibrium  at  Tju  (2600  K),  then  Xh.ju  =  0.21.  We  therefore  estimate 
the  value  of  Xh,o  to  be  0.21  x  0.011  =  2.4  x  10-3.  The  value  measured  by  Hsu  for 
these  conditions  is  2x  10~3  for  the  lowest-methane  case  reported  (0.5%)  [18].  Thus,  an 
estimate  based  on  simple  diffusion  theory  is  able  to  reproduce  the  experimental  result 
to  within  experimental  uncertainty.  (However,  an  uncertainty  of  at  least  a  factor  of 
two  should  be  assigned  to  the  estimated  value,  since  Xhju  is  probably  somewhat  less 
than  the  equilibrium  value  [16],  and  the  geometric  assumption  of  an  infinitely-long 
filament  and  a  cylinderical  substrate  will  introduce  some  error.) 

To  achieve  a  smaller  id  (and  therefore  a  larger  Xn,o )  in  a  filament  system  requires 
moving  the  filament  closer  to  the  substrate  or  using  a  larger  diameter  filament.  Taking 
these  ideas  to  the  logical  extreme,  a  conceivable  high-rate  diamond  growth  method 
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might  be  a  “sandwich”  geometry,  in  which  an  electrically-heated  refractory  metal 
mesh  is  separated  from  a  water-cooled  substrate  by  a  small  gap  (say  1  mm),  which 
contains  a  hydrogen /hydrocarbon  gas  mixture.  In  this  case,  id  is  approximately  the 
gap  width,  which  can  be  made  very  small,  and  film  uniformity  would  not  be  a  problem. 
While  the  electrical  power  requirements  would  be  large,  they  would  be  no  larger  than 
those  for  current  high-rate  methods  for  the  same  area. 

B.  Microwave  Plasma  Reactors 

Hsu  has  also  measured  the  H  mole  fraction  at  the  substrate  in  a  standard  low- 
Pe  microwave  plasma  system  at  20  Torr,  and  found  that  the  value  is  1.2  x  10~3 
—  close  to  that  for  the  filament  system  [19].  This  suggests  that  id  for  a  diffusion- 
dominated  microwave  system  is  also  several  centimeters.  Since  we  would  expect 
from  diffusion  theory  that  id  is  of  the  order  of  (but  somewhat  less  than)  the  plasma 
ball  dimensions,  this  estimate  appears  reasonable.  We  would  expect,  then,  that  H 
transport  is  usually  diffusion-limited  in  low-Pe  microwave  plasma  reactors,  in  the 
sense  that  Xh, o  Xh, plasma- 

Rau  and  Picht  [21]  have  reported  experimental  evidence  for  diffusion-limited  film 
growth  in  a  microwave  plasma  at  300  Torr,  by  growing  on  a  substrate  containing 
raised  flattened  cones.  Assuming  id  >  1  cm  in  this  experiment  (  a  conservative 
estimate),  we  estimate  using  the  above  procedures  that  Xh,o  I  Xh, plasma  <  5  x  10~3; 
therefore  H  transport  should  be  strongly  diffusion-limited.  Since  CH3  is  tightly- 
coupled  chemically  to  H  through  the  reaction  [3] 

CH<  +  H  ^  CH3  +  H2,  (15) 

the  CH3  distribution  will  follow  the  H  distribution.  Since  the  growth  rate  under  these 
conditions  is  expected  [1]  to  depend  on  both  [H]  and  [CH3],  the  spatial  dependence 
of  growth  rate  should  reflect  the  diffusion-limited  transport  of  H.  By  comparing  their 
experimental  results  with  simulations,  Rau  and  Picht  estimate  that  some  species 
which  affects  the  growth  rate  has  a  destruction  probability  on  the  surface  of  0.1. 
Although  they  assume  this  is  the  growth  species  (they  do  not  consider  H  transport), 
this  value  agrees  well  with  the  recombination  coefficient  of  H  on  diamond. 

These  examples  have  shown  that  simple  diffusion  theory  is  able  to  rationalize 
experimental  observations  in  both  filament  and  plasma  systems  when  Pe  <C  1.  The 
H  mole  fraction  at  the  substrate  can  be  estimated  with  reasonable  accuracy,  and 
the  observation  of  film  growth  features  associated  with  diffusion-limited  transport 
is  explained.  Diffusion  theory  also  indicates  the  modifications  needed  to  achieve 
higher  H  mole  fractions  at  the  surface  in  these  systems.  Higher  input  power  can  be 
used  to  fully  dissociate  hydrogen  near  the  filament  or  in  the  plasma;  beyond  this, 
it  is  necessary  to  decrease  id-  This  is  difficult  to  do  in  low-pressure,  low-Pe  plasma 
systems,  since  the  plasma  dimensions  are  governed  by  the  excitation  wavelength. 
However,  modifications  of  current  filament  systems  could  in  principle  achieve  a  small 
id  and  become  an  alternate  means  of  achieving  high  film  growth  rates. 
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IV.  CONVECTION-DOMINATED  METHODS 


Due  to  the  difficulties  of  achieving  a  small  diffusion  length  scale  at  low  Pe,  most 
high-rate  diamond  growth  techniques  make  use  of  convective  transport.  In  a  high-Pe 
(or  high-Re)  flow,  the  diffusion  length  scale  at  the  substrate  is  not  determined  by 
reactor  dimensions,  but  by  an  inherent  length  scale  of  the  flow  —  the  boundary  layer 
thickness. 

A  typical  configuration  for  high-rate  diamond  growth  is  shown  in  Fig.  1(a).  A 
high-velocity,  high-temperature  gas  jet  impinges  at  normal  incidence  on  a  cooled 
substrate  and  is  deflected  radially.  A  second  type  of  flow,  which  is  used  occasionally 
for  diamond  growth  [6,  7,  22,  23],  is  shown  in  Fig.  1(b).  Here  a  small  substrate  is 
immersed  in  a  larger,  uniform  gas  flow. 

The  flow  in  Fig.  1(a)  is  more  difficult  to  analyze  that  that  in  Fig.  1(b),  since  the 
effects  of  viscous  drag,  heat  conduction,  and  radial  species  diffusion  may  significantly 
broaden  and  cool  the  jet  after  it  leaves  the  torch  nozzle  before  impinging  on  the 
substrate.  For  this  reason,  we  will  consider  primarily  the  flow  in  Fig.  1(b),  although 
most  of  the  results  will  apply  qualitatively  to  the  flow  in  Fig.  1(a)  by  substituting 
dj  for  d4,  where  dj  is  an  approximate  jet  diameter  at  the  substrate  (or  the  deposit 
diameter).  We  will  assume  that  far  from  the  substrate  the  flow  is  uniform  and  is 
directed  downward  with  velocity  £/«,,  and  has  temperature  Tx  and  species  mole 
fractions  Xkt00- 

An  important  flow  parameter  is  the  approach  Mach  number  M^,,  defined  as  the 
ratio  of  £/<„  to  the  sound  speed  at  7! Both  subsonic  flows  (M,*,  <  1)  and  supersonic 
flows  (Mqo  >  1)  are  considered  here.  We  restrict  attention  to  the  continuum  flow 
regime,  where  the  mean  free  path  is  much  less  than  the  substrate  dimensions.  From 
Eq.  (10),  this  is  reasonable  for  pressures  above  a  few  Torr.  (While  most  convective 
diamond  growth  reactors  operate  in  this  flow  regime,  the  supersonic  arcjet  reactor  of 
boh  end  Cappelli  [24,25]  not,,  and  therefore  is  beyond  the  scope  of  this  analysis.) 

We  will  defer  the  issue  of  homogeneous  chemistry  to  section  IV  B,  where  re¬ 
combination  of  H  in  the  boundary  layer  is  examined,  and  begin  by  considering  the 
“frozen  flow”  limit  in  which  homogeneous  chemistry  is  neglected.  The  issue  of  homo¬ 
geneous  chemistry  outside  the  boundary  layer  (due,  for  example,  to  a  nonequilibrium 
freestream,  or  to  the  presence  of  a  shock  wave)  is  not  considered  here 

In  high  Reynolds  number  flows,  the  effects  of  diffusive  transport  (viscosity,  heat 
conduction,  and  molecular  species  diffusion)  are  significant  only  in  a  thin  region  near 
solid  surfaces  (the  boundary  layer)  and  in  the  wakes  behind  objects.  Within  the 
boundary  layer,  gradients  in  velocity,  temperature,  and  species  mole  fractions  are 
large:  the  radial  velocity  rapidly  decreases  from  its  value  outside  the  boundary  layer 
(given  by  potential  flow  theory)  to  zero  at  the  surface,  and  the  temperature  decreases 
from  Too  to  the  substrate  temperature.  (It  may  be  shown  from  kinetic  theory  that 
the  translational  gas  temperature  at  the  surface  must  equal  the  surface  temperature 
if  the  substrate  dimensions  are  much  larger  than  the  gas  mean  free  path  [26].  Wrhile 
temperature  discontinuities  have  been  measured  [27]  at  filaments  in  HFCVD  reactors, 
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this  is  possible  since  the  filament  diameter  is  similar  to  the  mean  free  path.)  If 
a  particular  species  reacts  on  the  surface  or  in  the  boundary  layer,  then  its  mole 
fraction  will  also  vary  sharply  through  the  boundary  layer. 

In  a  gas,  the  boundary  layers  for  velocity,  temperature,  and  species  are  all  similar 
in  thickness,  since  the  transport  properties  determining  each  (the  kinematic  viscosity 
v,  the  thermal  diffusivity  a,  and  the  diffusion  coefficients  £>*,  respectively)  are  all 
similar  in  magnitude.  The  boundary  layer  thickness  may  be  defined  in  various  ways. 
For  the  present  purposes,  the  most  useful  definition  is  the  conduction  thickness  (12). 
For  H  transport,  this  is  defined  as 


Sh 


Xh.qq  —  XH,o 

(dXH/dz)z- o 


(16) 


With  this  definition  for  the  boundary  layer  thickness,  it  is  clear  by  comparison  with 
Eq.  (7)  that  the  diffusion  length  scale  for  H  transport  is  simply  given  by  =  Sh.  A 
similar  thermal  boundary  layer  thickness  may  be  defined: 


f,  —  Too  ~  To 

T_  (dT/dz)z=0  ’ 


(17) 


The  geometric  interpretation  of  Sh  and  St  is  shown  in  Fig.  2. 

We  are  most  interested  in  the  boundary  layer  near  the  center  of  the  substrate  — 
i.e.,  near  the  stagnation  streamline.  Since  the  flow  decelerates  approaching  the  sub¬ 
strate,  the  velocity  is  always  small  in  this  region.  For  this  reason,  the  stagnation-point 
boundary  layer  is  always  laminar,  even  if  the  oncoming  flow  is  turbulent.  Neverthe¬ 
less,  turbulence  may  play  some  role,  since  it  is  known  that  freestream  turbulence  can 
enhance  heat  and  mass  transfer  through  the  laminar  stagnation  point  boundary  layer 
[28,  29).  The  possible  effects  of  freestream  turbulence  will  not  be  considered  here. 

The  boundary  layer  thickness  at  the  stagnation  point  is  determined  by  the  trans¬ 
port  properties  and  a  single  parameter  a  which  characterizes  the  velocity  field  outside 
the  boundary  layer.  The  a  parameter,  known  either  as  the  stagnation-point  velocity 
gradient  parameter  or  the  strain  rate,  is  defined  by 


a  =  dv/dr. 


(18) 


where  v  is  the  radial  (tangential)  velocity  just  outside  the  boundary  layer.  For  radii 
not  too  large,  a  is  independent  of  r  [26],  and  thus  the  boundary  layer  thickness  is 
constant  over  the  central  portion  of  the  substrate.  It  will  be  shown  in  Section  IV  R 
that  a-1  is  the  characteristic  time  for  species  to  diffuse  through  the  boundary  layer, 
and  therefore  a  may  be  thought  of  also  as  an  inverse  boundary  layer  diffusion  time. 
From  dimensional  considerations,  the  a  parameter  must  scale  with  the  ratio  U^/d,: 

a  =  C~.  (19) 
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The  constant  C  depends  <«n  the  geometry  of  the  substrate  (i.e.,  disk,  sphere,  cylinder, 
etc.),  as  well  as  the  .•’-proach  Mach  number  of  the  flow  Moo.  Values  for  some  simple 
geometries  are  given  in  Table  II. 

If  the  temperature  difference  between  the  freestream  and  the  surface  is  small,  so 
that  the  transport  properties  and  density  are  constant  through  the  boundary  layer, 
the  .  simple  analytical  expressions  for  the  boundary  layer  thickness  may  be  obtained 
[12].  At  the  stagnation  point, 


and 


(20) 


(21) 


assuming  ajv  and  D^jv  are  not  too  far  from  1.0. 

For  a  dilute  hydrocarbon  mixture  in  H?  at  1200  K,  v  ~  8  x  10 3/p  cm2/s,  and 
o  w  1.2  x  104/p  cm2/s  ( p  in  Torr)  [8].  Therefore, 


and 


f  140 
bj  22  — —  cm 
y/ap 


(22) 


r  160 
oh  ~  ~r=  cm- 

s/ap 


(23) 


with  a  in  s-1  and  p  in  Torr. 

Although  Eqs.  (22)  and  (23)  are  derived  for  the  case  of  constant  properties  (small 
temperature  changes),  comparison  with  numerical  simulation  results  shows  that  they 
are  in  fact  reasonably  accurate  for  typical  boundary  layers  in  diamond  CVD  reactors, 
with  temperature  changes  of  up  to  4000  K  through  the  boundary  layer.  For  the  sim¬ 
ulations,  the  gas  was  assumed  to  be  a  1%  CH4  in  H2  mixture  in  chemical  equilibrium 
outside  the  boundary  layer,  and  chemistry  in  the  boundary  layer  was  neglected.  Both 
multicomponent  mass  diffusion  and  thermal  diffusion  were  included.  The  boundary 
layer  thicknesses  from  the  simulations  agreed  with  those  predicted  by  Eqs.  (22)  and 
(23)  to  within  ±20%  for  freestream  temperatures  from  3000  K  to  5000  K. 

Using  Eq.  (23)  for  in  Eq.  (11)  results  in 


AVo 

A//.00 


1  ±  11007// 


p  (Torr) 
a  (s-1) 


(24) 


Equation  (24)  shows  that  H  transport  through  the  boundary  layer  is  determined 
by  the  ratio  pja.  For  optimal  H  transport,  the  flow  velocity  (or  Mach  number) 
should  be  chosen  to  give  an  a  value  such  that  the  second  term  in  Eq.  (24)  is  of 
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order  1;  for  smaller  a ,  H  transport  will  be  diffusion-limited  and  therefore  H  will 
be  significantly  depleted  at  the  substrate,  while  larger  a  values  do  not  substantially 
increase  Xh.o-  Therefore,  since  7#  %  0.1,  p/a  should  be  of  order  10~4  Torr-s  or  less 
to  avoid  significant  diffusion  limitations  to  H  transport. 

Simulations  indicate  that  typical  oxyacetylene  torch  flames  at  1  atm  operate  with 
a  in  the  range  103  to  104  s-1  [5].  Therefore  p/a  as  0.08  to  0.8  Torr-s,  which  is  well 
into  the  diffusion-limited  regime.  On  the  other  hand,  some  plasma  jet  reactors  may 
achieve  a  values  as  large  as  10s  to  106  s-1  for  small  substrates  (several  mm)  (5).  In 
this  case,  p/a  is  between  2  x  10~4  and  2  x  10-3  Torr-s,  near  the  edge  of  the  diffusion- 
limited  regime.  For  larger  substrates,  however,  it  is  likely  that  here  too  H  transport 
is  nearly  completely  diffusion  limited. 

In  the  diffusion-limited  regime,  Eq.  (24)  reduces  to 


XHt0  9.1  x  10 


-4 


a  (*) 


Xh, 00  7 H  p  (Torr) 

Since  [H]0  =  XHi0(p/ RT0),  for  T0  =  1200  K 

1  2  x  10_n 

[H]0  = - A 'h, coJap  moles/cm3, 

‘7  H 


(25) 


(26) 


with  p  in  Torr  and  a  in  s_1.  This  expression  shows  that  if  Xh,oo  is  constant  and 
boundary-layer  chemistry  may  be  neglected,  then  the  H  concentration  at  the  substrate 
in  the  diffusion-limited  regime  is  proportional  to  y/ap.  Therefore,  increasing  either  a 
or  p  will  result  in  an  increase  in  [H]o-  However,  as  will  be  discussed  in  the  next  section, 
a  cannot  be  increased  indefinitely.  Also,  at  sufficiently  high  pressures  H  will  begin 
recombining  homogeneously  in  the  boundary  layer,  which  will  set  an  upper  limit  on 
the  desirable  operating  pressure.  In  the  sections  below,  we  turn  to  the  question  of  the 
“optimal”  conditions  which  maximize  the  H  concentration  at  the  surface.  It  will  be 
shown  that  there  exists  both  an  optimal  flow  Mach  number  and  an  optimal  pressure. 


A.  The  Optimal  Mach  Number  for  H  Transport 

Since  a  =  CU^/d,,  to  increase  a  for  a  specified  substrate  diameter  requires  in¬ 
creasing  the  flow  velocity  U0 0.  However,  there  are  limits  on  the  velocity  to  which  a 
gas  may  be  accelerated,  since  the  flow  acceleration  is  provided  by  expanding  the  gas 
through  a  nozzle.  The  largest  obtainable  velocity  is  of  the  order  of  the  sonic  velocity 
at  the  temperature  upstream  of  the  nozzle.  Therefore,  if  this  temperature  is  fixed  (for 
example,  by  materials  constraints),  then  there  is  a  maximum  attainable  flow  velocity, 
and  therefore  a  maximum  dsa. 

To  get  an  estimate  of  amax,  let  us  consider  a  gas  with  constant  specific  heat, 
which  is  at  temperature  T°  and  pressure  p°  upstream  of  the  nozzle,  and  is  expanded 
reversibly  and  adiabatically  to  a  Mach  number  Moo.  (1°  a  real  nozzle,  heat  losses  and 
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other  irreversibilities  are  significant.  In  this  case,  T°  and  p°  should  be  interpreted 
as  the  stagnation  temperature  and  pressure,  respectively,  downstream  of  the  nozzle, 
which  may  be  significantly  less  than  the  upstream  static  temperature  and  pressure.) 
In  this  case  [30], 


U oo  —  Moo 


pry/2 


(27) 


where  k  =  Cp/c*,  and  m  is  the  molecular  weight  of  the  gas.  Note  that  i'oc  asymptot 
ically  approaches  a  finite  value  as  Moo  — *  oo. 

Since  a  —  C(M00)U00/d3 ,  we  may  substitute  from  Eq.  (27)  for  resulting  in 

1  f  RT°\  1/2 
°  “  d,  {  m  ) 


/4(Moo,fc) 


(28) 


where 


>i(Moo,  &)  =  CMc 


V  +  (¥)M2oo' 


(‘29) 


The  function  y4(Moo,fc)  is  plotted  in  Fig.  3(a)  for  the  two  k  values  1.3  and  1.7,  which 
bracket  those  expected  for  a  mixture  of  H  and  H2.  These  calculations  were  done 
for  a  spherical  substrate,  since  data  are  available  on  the  dependence  of  C  on  Mach 
number  for  this  geometry.  ( C  decreases  from  3  at  low  Mach  number  to  1.14  at  high 
Mach  number  [26].  These  data  take  into  account  the  presence  of  the  bow  shock  for 

Moo  >  1.) 

The  results  in  Fig.  3(a)  show  that  ^(M^,*:)  attains  a  limiting  value  of  approx¬ 
imately  3  at  high  Mach  numbers.  For  a  flat  disk,  the  dependence  of  C  on  Mach 
number  is  not  as  well  known;  however,  C  is  estimated  to  be  about  a  factor  of  three 
lower  than  for  a  sphere  at  high  Mach  numbers  [26].  The  maximum  value  of  a  for  a 
disk  is  then 


_  j_  (  rt°v  12 

d3\m) 


(30) 


From  this  relation  it  is  clear  that  high  T°  and  low  gas  molecular  weight  result  in  the 
largest  possible  velocity-gradient  parameter.  If  we  take,  as  extreme  values,  T°  —  104 
K  and  m  =  1,  then  amaxd3  «  106  cm/s.  Therefore,  if  d,  =  10  cm,  then  amaz  105  s-1. 
From  Eq.  (24),  for  this  substrate  diameter  H  transport  will  be  necessarily  diffusion- 
limited  for  pressures  greater  than  approximately  10  Torr.  Since  most  high-rate  growth 
methods  operate  at  a  pressure  significantly  higher  than  this,  this  calculation  illustrates 
that  in  general,  large-area  growth  processes  always  operate  in  the  diffusion-limited 
regime. 

The  effect  of  molecular  weight  on  amax  is  worth  noting,  since  some  plasma  methods 
make  use  of  an  argon /hydrogen  mixture.  If  we  substitute  Eq.  (21)  into  Eq.  (9)  and 
use  Eq.  (30),  we  find  approximately  that  in  the  diffusion  limit, 
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(31) 


If  the  mean  molecular  weight  of  the  gas  is  increased  by  argon  addition,  both  Dh 
and  amai  decrease  (since  the  sound  speed  is  lowered),  and  therefore  the  maximum  H 
concentration  achievable  at  the  substrate  is  lower  (in  addition  to  the  effects  of  argon 
dilution  on  Xh,oo)- 

In  the  absence  of  homogeneous  chemistry,  the  maximum  H  concentration  at  the 
surface  will  be  achieved  when  the  product  ap  is  maximized,  where  p  is  the  pressure 
at  the  substrate.  Therefore,  the  variation  of  p  with  Mach  number  must  al  be 
considered.  For  a  subsonic  flow,  pressure  recovery  is  nearly  complete,  and  therefore 
p  «  p°.  For  a  supersonic  flow,  however,  a  shock  wave  forms  above  the  substrate.  Due 
to  the  loss  of  stagnation  pressure  across  the  shock  wave,  p  at  the  substrate  will  be 
less  than  p°  if  >  1.  Therefore,  for  a  given  upstream  stagnation  pressure  p°, 

[H]o  cx  \/a(p/p0)  (32) 

where  p/p°  is  the  stagnation  pressure  ratio  across  the  shock.  In  Fig.  3(b),  the  de¬ 
pendence  of  [H]o  on  Moo  given  by  Eq.  (32)  is  shown.  For  these  calculations,  a  is 
taken  from  Fig.  3(a),  and  p/p°  is  calculated  from  standard  normal  shock  relations 
[30],  These  results  show  that  [H]0  is  maximized  at  a  Mach  number  between  1  and  2 
for  the  frozen  flow  case  (no  chemistry). 

In  reality,  the  frozen  flow  assumption  may  not  hold,  particularly  if  a  shock  wave  is 
present.  At  the  high  temperature  and  pressure  behind  the  shock,  the  H  mole  fraction 
may  equilibrate  rapidly;  in  this  case,  the  value  of  Xh,oo  (behind  the  shock  but  outside 
the  boundary  layer)  may  be  dependent  on  Mach  number.  Therefore  the  Mach  number 
at  which  [H]0  is  maximized  may  differ  somewhat  from  that  for  the  frozen  flow  case. 
To  better  determine  the  optimal  Mach  number  would  require  a  full  reacting  flow 
simulation  of  the  shock  layer,  which  is  beyond  the  scope  of  the  present  work. 


B.  Homogeneous  Recombination 


In  addition  to  recombination  of  atomic  hydrogen  on  the  diamond  surface,  H  may 
recombine  in  the  gas-phase  at  high  pressures.  This  too  will  affect  H  at  the  surface, 
and  will  set  an  upper  limit  on  the  pressure  at  which  it  is  desirable  to  operate. 

To  examine  the  effects  of  homogeneous  recombination,  let  us  consider  a  flow 
which,  outside  the  boundary  layer,  is  in  chemical  equilibrium  at  T^.  In  this  case,  the 
H  concentration  outside  the  boundary  layer  is  maintained  by  the  temperature,  and 
recombination  only  occurs  within  the  boundary  layer. 

The  characteristic  time  for  an  H  atom  at  the  edge  of  the  boundary  layer  to  diffuse 
to  the  substrate  is  approximately 


Td  2S 


SJL. 

Dh ' 


(33) 
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Using  Eq.  (21)  for  6h  results  in 


Tjss  a  1 . 


(34) 


Therefore,  the  a  parameter  may  be  thought  of  as  an  inverse  boundary  layer  diffusion 
time,  in  addition  to  its  fluid-mechanical  interpretation  as  a  velocity  gradient. 

We  may  also  define  an  effective  first-order  recombination  rate  coefficient  an  as 


m 

dt 


(35) 


Clearly,  homogeneous  recombination  in  the  boundary  layer  will  be  important  when¬ 
ever  Or  >  a.  The  nondimensional  ratio  a/*/a  defines  the  Damkohler  number  Da^  for 
this  problem.  Recombination  is  therefore  significant  when  Da^  >  1 . 

In  a  pure  hydrogen  gas,  H  recombines  through 

H  +  H  +  H2  -►  2H2.  (36) 

In  this  case, 

a*  =  *a«[H](HaJ,  (37) 

where  k 36  is  the  rate  constant  for  this  reaction,  equal  to  9.7  x  10-6T~06  cm6/(mo]2  s) 
[31].  At  constant  mole  fraction,  Eq.  (37)  shows  that  or  oc  p2.  Therefore,  for  a  given 
degree  of  dissociation,  the  condition  Daj  =  1  will  be  attained  when  the  parameter 
group  p/a1/2  exceeds  some  critical  value. 

With  a  small  amount  of  added  hydrocarbon,  a  second  path  competes  with  reaction 
(36),  and  in  many  cases  dominates  the  recombination  rate.  This  path  is  due  to  the 
two  reactions 


CH3  +  H  +  M  -  CH4  +  M  (38) 

and 

CH4  +  H  ^  CH3  +  H2.  (39) 

The  rate-limiting  step  is  reaction  (38),  which  is  in  the  low-pressure  (termolecular) 
limit  even  for  pressures  greater  than  1  atm  [32],  and  is  effectively  irreversible.  As 
soon  as  a  methyl  is  consumed  by  reaction  (38),  an  H  atom  is  abstracted  from  the 
resulting  methane  molecule  by  (39),  maintaining  a  steady-state  methyl  concentration, 
but  resulting  in  a  net  loss  of  H  atoms.  This  path  is  fast,  since  reaction  (38)  has  no 
activation  barrier  [32]  and  the  barrier  for  reaction  (39)  is  less  than  9  kcal/mole  [31]. 

Since  the  CH3  and  CH4  concentrations  in  the  boundary  layer  are  strongly  affected 
by  boundary-layer  chemistry,  it  is  difficult  to  develop  criteria  for  the  significance  of 
homogeneous  H  recombination  based  on  simple  arguments  when  hydrocarbons  are 
present.  In  this  case,  it  is  preferable  to  examine  the  results  of  numerical  simulation 
of  the  boundary  layer. 
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As  an  example,  calculated  species  profiles  are  shown  in  Fig.  4  for  a  set  of  conditions 
where  H  recombination  is  significant.  For  these  conditions,  acetylene  is  the  most 
abundant  hydrocarbon,  and  is  hardly  affected  by  boundary- layer  chemistry.  Methane, 
on  the  other  hand,  increases  by  two  orders  of  magnitude  within  a  few  mm  of  the 
substrate.  The  H  profile  is  seen  to  be  strongly  affected  by  boundary-layer  chemistry. 

The  contributions  of  several  reactions  to  the  homogeneous  recombination  rate  for 
the  conditions  of  Fig.  4  are  shown  in  Fig.  5.  Reaction  (36)  contributes  24%  to  the 
total  recombination  rate  integrated  through  the  boundary  layer  for  these  conditions, 
and  reactions  (38)  and  (39)  each  contribute  34%.  Therefore,  these  three  reactions 
together  account  for  92%  of  the  total  homogeneous  recombination  rate.  Although 
acetylene  is  the  most  abundant  hydrocarbon,  the  reaction 


H  +  C2H2  +  M  -»  C2H3  +  M 


(40) 


only  contributes  3%  to  the  H  recombination  rate  under  these  conditions. 

To  quantify  the  effect  of  homogeneous  recombination,  a  correction  factor  Fr  may 
be  defined  as 

r  _ 

'  R  "  (A 


where  (X//,0)nc  is  the  H  mole  fraction  at  the  substrate  neglecting  boundary  layer 
chemistry  (but  still  including  heterogeneous  recombination  on  the  substrate).  The  H 
concentration  at  the  substrate  is  then  determined  by  multiplying  Eq.  (26)  by  Fr. 

Calculations  of  boundary  layer  profiles  were  carried  out  with  and  without 
boundary-layer  chemistry  for  a  range  of  a  and  p  values,  and  the  correction  factor 
Fr  determined.  The  gas  was  taken  to  be  a  mixture  of  1%  CH4  in  H2,  and  the  gas 
outside  the  boundary  layer  was  assumed  to  be  in  chemical  equilibrium  at  Too  =  3000 
K.  The  factor  Fr  is  shown  in  Fig.  6  plotted  against  p/a1/2  for  the  three  pressures  100, 
1000,  and  10000  Torr.  It  is  seen  that  the  dependence  of  Fr  on  p/a1/2  is  approximately 
the  same  for  all  three  pressures,  and  therefore  pi  a}12  is  a  valid  approximate  scaling 
parameter  for  homogeneous  recombination.  Homogeneous  recombination  is  seen  to 
become  significant  for  p/a1/2  greater  than  a  few  Torr-s1/2. 

Since  Fr  falls  off  approximately  linearly  with  pressure  in  the  regime  where  re¬ 
combination  is  significant,  while  in  the  absence  of  recombination  [H]o  oc  p1/2,  for  a 
specified  a  and  XfftO0  there  will  be  some  pressure  at  which  [H]o  is  maximized,  above 
which  it  falls  off  as  p-1/2.  If  Xri00  also  varies  with  pressure  (as  it  may  if  the  gas 
outside  the  boundary  layer  is  in  chemical  equilibrium),  then  the  scaling  with  p  will 
be  modified,  but  a  maximum  at  a  particular  pressure  is  still  expected. 

This  behavior  is  shown  in  Fig.  7,  in  which  simulation  results  for  [H]o  are  plotted 
vs.  p  at  fixed  a  and  Tx,  assuming  chemical  equilibrium  outside  the  boundary  layer. 
For  Too  =  3000  K,  [H]o  is  nearly  independent  of  p  over  a  wide  range  if  boundary 
layer  chemistry  is  neglected,  due  to  the  decline  of  the  equilibrium  value  of  Xr.oc  with 
increasing  p.  With  chemistry  included,  [H]o  peaks  and  then  decreases  at  high  p.  The 
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pressure  at  which  [H]o  is  maximal  is  seen  to  be  just  where  homogeneous  recombination 
begins  to  become  significant,  given  approximately  by  p/a1(,a  =  2  Torr-s1/2. 

For  Too  =  5000  K,  the  freestream  is  nearly  fully  dissociated  up  to  103  Torr,  and  the 
pressure  dependence  is  approximately  as  p1^2  at  low  p.  For  this  more  activated  gas, 
[Hjo  is  maximized  at  p/a1^2  «  5  Torr-s1/2.  For  a  =  10s  s-1,  the  maximum  value  of  [H]0 
is  seen  to  be  approximately  10~6  moles/cm3.  Since  this  a  value  is  the  largest  achiev¬ 
able  for  a  10  cm  diameter  substrate,  and  at  5000  K  the  freestream  is  fully  dissociated, 
this  H  concentration  is  an  absolute  upper  bound  on  the  H  concentration  achievable 
by  means  of  convective  transport  at  the  surface  of  a  10  cm  diameter  substrate.  The 
pressure  at  which  this  maximum  occurs  is  approximately  2  atmospheres. 

C.  Substrate  Heat  Flux 

The  conditions  which  lead  to  high  H  concentrations  at  the  surface  also  result  in 
substantial  heat  fluxes  to  the  substrate,  which  may  in  fact  be  a  limiting  consideration 
for  process  design.  Here  we  develop  an  expression  for  the  substrate  heat  flux  for  the 
case  where  boundary  layer  H  recombination  can  be  neglected  (p/a1^2  <  1  Torr-s1/2). 

Substrate  heating  results  primarily  from  two  sources:  thermal  conduction,  and 
heating  due  to  recombination  of  atomic  hydrogen.  The  heat  flux  to  the  substrate  may 
be  written  as 

q,  =  k[in)  <42) 

where  k  is  the  thermal  conductivity  of  the  gas  at  the  substrate,  A  Hr  is  the  heat  of 
recombination  of  atomic  hydrogen  (223  kJ/mole  at  1200  K)  and  Rh  is  the  atomic 
hydrogen  recombination  rate  (moles  of  H/cm2/s).  From  Eq.  (8), 

Rh  =  no~r^-  {Xh,oo  —  Xn,o)  ■  (43) 

*>H 

Using  the  definition  of  Eq.  (42)  may  be  rewritten  as 

(44) 

where  AT  =  Too  —  To  and  AX#  =  Xh>00  —  Xhi0.  In  writing  Eq.  (44),  we  have  also 
made  use  of  the  relation  k  =  anocp,  where  cp  is  the  molar  heat  capacity  at  constant 
pressure  at  the  substrate.  Since  the  ratio  6t/6h  is  simply  ( a/ Dh  )0,4  [Eqns.  (20)  and 
(21)],  we  may  simplify  this  to 

=  ^  cpAT  +  A//fl(^)°'6AXH  .  (45) 

Substituting  numerical  values  we  have 
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qa  =  3.5  x  KT5^]  AT +  8600 AX/,],  (46) 

where  qa  is  in  W/cm2,  a  is  in  s-1,  p  is  in  Torr,  and  AT  is  in  K.  The  heat  flux  given  by 
this  expression  agrees  with  simulation  results  to  within  about  10%.  From  Eq.  (46), 
it  is  clear  that  heat  conduction  and  recombination  heating  will  contribute  similar 
amounts  to  the  total  substrate  heat  flux,  since  AT  may  be  of  order  2000  K,  and 
AXh  <  1. 

Since  both  the  heat  flux  and  the  H  concentration  at  the  surface  depend  on  y/ap, 
a  high  H  concentration  at  the  surface  will  always  be  accompanied  by  a  large  heat  flux 
to  the  substrate.  Some  typical  values  illustrate  the  magnitudes  involved.  If  we  take 
a  —  10s  s-1,  p  —  1  atm,  and  X//)00  —  0.5,  then  q,  =  1300  W/cm2. 

For  comparison,  the  heat  flux  to  the  substrate  during  typical  oxyacetylene  torch 
growth  is  roughly  200  W/cm2,  while  the  heat  flux  to  the  water-cooled  cathode  in  an 
arcjet  is  of  order  2000  W/cm2  [33].  For  a  more  extreme  case,  the  heat  flux  at  the 
throat  of  the  nozzle  of  the  Space  Shuttle  main  engine  is  of  order  13,000  W/cm2  [34]. 
In  this  case,  cooling  is  accomplished  using  integral  cooling  passages  in  the  nozzle  wall 
which  carry  pressurized  hydrogen. 

These  heat  flux  estimates  show  that  an  important  component  of  reactor  design 
for  high  growth  rates  is  the  substrate  thermal  design.  To  handle  these  heat  fluxes  and 
achieve  reproducible  substrate  temperatures,  it  is  critical  to  eliminate  thermal  contact 
resistance  and  low-conductivity  materials.  Therefore,  metal  substrates  (for  example, 
molybdenum)  with  integral  cooling  passages  will  be  required  for  high-rate,  large-area 
diamond  growth.  Since  the  thermal  conductivity  of  molybdenum  at  1200  K  is  about  1 
W/cm/K  [35],  the  temperature  drop  through  a  1  mm  thick  sheet  would  be  130  K  for  a 
heat  flux  of  1300  W/cm2.  Therefore,  molybdenum  would  be  an  acceptable  substrate 
material.  Pressurized  boiling  water  would  probably  be  an  acceptable  coolant.  With 
such  a  substrate  design,  heat  fluxes  of  order  1000-2000  W/cm2  should  not  present 
major  problems. 


V.  PROCESS  SCALING  RELATIONS 

The  results  in  the  previous  sections  may  be  combined  to  yield  relations  relating 
growth  rate  and  process  efficiency  to  controllable  process  parameters.  The  most 
important  relation  for  process  scaling  is  Eq.  (1),  which  may  be  rewritten  as 

(47) 

For  a  specified  defect  density  (film  quality),  the  sustainable  growth  rate  is  predicted 
to  scale  quadratically  with  [H]o-  Of  course,  the  model  which  led  to  this  expression  was 
highly  simplified  [1].  For  this  reason,  this  quadratic  dependence  on  [H]0  is  uncertain. 
Experimental  data  on  the  relation  of  G ,  X*/,  and  [H]o  are  needed;  in  the  absence  of 
such  data,  here  we  explore  the  scaling  implications  of  Eq.  (47). 

In  the  regime  where  H  transport  is  diffusion  limited  but  boundary  layer  recom¬ 
bination  is  insignificant,  [H]0  oc  •v/apX/,)0o  [Eq.  (26)].  Using  this  result  in  Eq.  (47) 
results  in 
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G  oc  (ap)XdeJX]fi00.  (48) 

Therefore,  to  maximize  G  for  a  specified  Xdtj ,  the  velocity  gradient  parameter  should 
be  maximized  (sonic  flow),  and  the  pressure  at  the  substrate  should  be  set  at  that 
value  where  homogeneous  recombination  just  begins  to  become  important,  as  dis¬ 
cussed  above.  Even  more  important,  the  gas  should  be  highly  activated,  since  the 
growth  rate  for  a  specified  Xdej  is  quadratic  in  the  degree  of  gas  activation  (Xh,oo)- 

Note  that  since  a  oc  l/d»,  the  growth  rate  at  constant  defect  density  is  inversely 
proportional  to  substrate  diameter.  Therefore,  growing  over  a  larger  substrate  diame¬ 
ter  at  the  same  rate  and  quality  as  achieved  for  small  d„  requires  either  increasing  the 
flow  velocity  (if  the  flow  is  subsonic)  to  maintain  the  same  a,  increasing  the  pressure, 
or  achieving  higher  gas  activation. 

This  also  suggests  that  if  it  is  not  necessary  to  grow  a  large-area  film,  then  it  is 
preferable  to  mount  several  small-area  substrates  in  the  gas  flow  (such  that  gas  can 
flow  between  the  substrates)  rather  than  a  single  substrate  of  equivalent  area.  In 
this  case,  each  small  substrate  will  have  a  larger  a  parameter  (and  therefore  thinner 
boundary  layer)  than  would  the  single  large  substrate. 

The  process  efficiency  may  be  more  important  than  the  absolute  growth  rate  in 
many  cases.  Several  different  measures  of  efficiency  may  be  defined  -  for  example, 
the  ratio  of  mass  deposition  rate  of  diamond  to  the  gas  flow  rate  (r/g),  or  the  ratio 
of  mass  deposition  rate  to  the  electrical  power  required  (r)e).  The  mass  flow  rate  per 
unit  area  may  be  expressed  in  terms  of  the  upstream  stagnation  temperature  T°  and 
pressure  p° ,  Mach  number  M^,  and  specific  heat  ratio  k  [33]. 

Assuming  that  Eq.  (48)  is  valid,  and  using  Eqs.  (28)  and  (29)  for  a,  we  find  that 


1 h  K 


rpQ 


ds 


(49) 


where,  as  above,  p/p°  is  the  stagnation-pressure  ratio  across  the  shock  if  >  1 
and  is  1.0  otherwise.  This  function  is  shown  in  Fig.  8,  normalized  to  the  value  at  low 
Mach  number.  The  efficiency  t?s  is  higher  for  supersonic  conditions,  since  although 
the  growth  rate  is  low  here,  the  gas  density  (and  therefore  the  gas  flow  rate  per  unit 
area)  is  even  lower.  However,  the  increased  efficiency  (a  factor  of  2)  may  not  be 
enough  to  justify  operation  under  supersonic  conditions. 

Equation  (49)  shows  that  is  strongly  dependent  on  Xu, <*>;  therefore,  the  disso¬ 
ciation  degree  should  be  maximized  for  high  process  efficiency.  The  process  efficiency 
is  independent  of  pressure,  if  H  transport  is  diffusion-limited,  since  both  the  growth 
rate  and  the  gas  flow  rate  are  linear  in  pressure.  The  process  efficiency  is  also  in¬ 
versely  proportional  to  d,  -  again,  this  suggests  that  there  is  a  price  to  be  paid  for 
large  substrate  diameter. 

Similar  scaling  relations  would  hold  for  t]e.  If  it  is  assumed  that  power  consump¬ 
tion  is  due  largely  to  the  energy  required  to  dissociate  hydrogen,  then  T]e  oc  J?g/Aj/i00. 
Therefore,  the  same  scaling  with  pressure,  substrate  diameter,  and  Mach  number  is 
obtained,  but  qe  is  now  only  linearly  proportional  to  X//iOC. 
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It  should  be  noted  that  these  conclusions  stem  from  Eq.  (47)  -  in  particular,  from 
the  result  that  the  growth  rate  at  constant  defect  density  scales  faster  than  linear 
with  [H]o-  If,  instead,  G  scaled  only  linearly  with  [H]0,  then  these  conclusions  would 
be  modified  significantly.  In  particular,  in  this  case  the  greatest  efficiency  would  be 
achieved  at  the  lowest  pressures.  Clearly,  experimental  data  on  the  scaling  of  growth 
rate  at  constant  defect  density  with  [H]o  are  needed. 

The  efficiencies  discussed  here  are  “single-pass”  efficiencies,  which  do  not  account 
for  possible  gas  or  energy  cycling.  Employing  cycling  is  the  most  significant  way  to 
increase  the  overall  process  efficiency.  Hydrogen  may  in  principle  be  indefinitely  recy¬ 
cled,  and  the  hydrocarbon  feedstock  need  only  be  added  to  make  up  that  lost  due  to 
film  growth.  Also,  the  heat  removed  from  the  substrate  and  other  reactor  components 
(cathodes,  walls)  can  be  used  to  preheat  the  feedstock  gas,  or  to  cogenerate  electrical 
power. 


VI.  SUMMARY  AND  CONCLUSIONS 

In  this  paper,  scaling  relations  have  been  developed  for  atomic  hydrogen  transport 
in  diamond  chemical  vapor  deposition  reactors.  It  was  shown  that  simple  ideas  based 
on  diffusion  theory  suffice  to  estimate  the  H  concentration  at  the  substrate  in  low- 
Pe  hot-filament  and  microwave  plasma  systems.  For  convection-dominated,  high-Pe 
reactors,  H  transport  is  dominated  by  the  characteristics  of  the  boundary  layer.  It  was 
shown  that  most  reactors  operate  in  a  regime  in  which  H-transport  is  diffusion-limited, 
in  which  case  the  H  concentration  at  the  substrate  scales  with  ^/dp,  if  homogeneous 
recombination  may  be  neglected.  Loss  of  H  through  homogeneous  recombination  in 
the  boundary  layer  was  shown  to  occur  primarily  through  two  chemical  pathways 
(direct  recombination,  and  through  CH3  +  H  +  CH4  +  M).  Since  these  reactions 
are  termolecular,  the  effects  of  homogeneous  recombination  were  shown  to  depend  on 
p/a1/2,  and  to  become  significant  when  this  parameter  exceeds  a  few  Torr-s1/2. 

Finally,  the  optimal  conditions  which  maximize  H  transport  where  determined 
for  high-Pe  reactors.  It  was  shown  that  approximately  sonic  flow  at  a  pressure  of 
roughly  1  or  2  atmospheres  is  optimal  for  H  transport,  and  therefore  should  allow  the 
highest  film  growth  rates.  Achieving  these  conditions  is  difficult  in  current  laboratory- 
scale  CVD  reactors,  due  to  materials  constraints  and  heat  losses.  It  may  be  possible 
to  achieve  them  in  large-scale  (megawatt)  arcjets,  or  possibly  in  inductively-coupled 
plasma  torches  or  microwave  torches.  Also,  chemical  rockets  might  provide  a  similar 
hot,  high-speed  gas  flow,  with  additional  benefits  due  to  residual  oxygen. 

There  are  of  course  many  issues  not  considered  here  which  will  affect  process 
scaling  and  optimization.  For  example,  in  plasma  reactors  lower  pressures  may  be 
necessary  for  discharge  stability  and  flow  uniformity  over  reasonable  areas.  Also, 
there  are  economic  trade-offs  which  must  be  assessed  in  comparing  different  reactor 
designs  to  determine  conditions  which  minimize  cost.  For  example,  achieving  the 
highest  possible  growth  rate  will  not  lead  to  the  minimum  cost  if  the  reactor  must 
be  significantly  more  complex  than  a  simpler  design  which  has  a  somewhat  lower 
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growth  rate.  Nevertheless,  the  relations  presented  here  provide  a  starting  point  both 
for  optimizing  current  reactors  and  designing  larger-scale,  next-generation  reactors. 
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TABLES 


TABLE  I.  Diffusion  length  scales  in  the  low-Pe  limit,  assuming  constant  diffusion  coef¬ 
ficient  and  no  gas- phase  chemistry. 


Geometry 

f-d,\ 

U,2 

Infinite  parallel  planes  separated  by  L 

L 

L 

Infinitely-long  concentric  cylinders,  R\  <  R2 

R1\n(R2/Rl) 

R2  ln(R2/ R\ ) 

Concentric  spheres,  R\  <  R2 

R1(l-R1/R2) 

R2{Ri/Ri-1) 

TABLE  II.  The  constant  C  in  Eq.  (19)  for  flov'  rormal  to  spheres,  cylinders,  and  disks 
of  diameter  Da  in  the  limits  of  low  and  high  approach  Mach  number.  [26,36] 

C 

Substrate  geometry  M, 

30  — *  0 

Moo  -*•  oo 

Sphere 

3 

1.14 

Cylinder 

4 

1.14 

Disk 

4/?r 

~  0.3 
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FIGURES 


FIG.  1.  Typical  configurations  for  high-rate  diamond  growth,  (a)  a  small  high-velocity 
jet  impinging  on  a  larger  substrate;  (b)  a  small  substrate  immersed  in  a  uniform  gas  flow. 


FIG.  2. 

6f. 


Geometric  interpretation  of  the  boundary  layer  conduction  thicknesses  6#  and 


FIG.  3.  (a)  Dependence  of  the  velocity  gradient  parameter  for  flow  past  a  sphere  on 

Mach  number;  (b)  Relative  H  concentration  at  the  substrate  vs.  Mach  number.  The  solid 
curves  are  for  specific  heat  ratio  k  —  1.3,  and  the  dashed  curves  for  k  =  1.7. 


FIG.  4.  Selected  species  profiles  calculated  for  a  1%  methane/hydrogen  mixture,  as¬ 
suming  chemical  equilibrium  outside  the  boundary  layer,  with  parameter  values  T, oo  —  3000 
K,  a  =  1000  s-1,  p  =  1000  Torr,  To  =  1200  K,  and  7//  =  0.1.  The  dashed  line  shows  the  H 
profile  resulting  if  boundary  layer  chemistry  is  neglected. 


FIG.  5.  Contribution  factors  for  several  reactions  to  H  recombination,  for  the  condi¬ 
tions  of  Figure  4. 


FIG.  6.  Ratio  of  H  mole  fraction  at  the  substrate  including  boundary  layer  chemistry 
to  that  calculated  neglecting  boundary  layer  chemistry  vs.  p/a U2.  The  conditions  other 
than  p  and  a  are  the  same  as  those  for  Figure  4. 


FIG.  7.  H  concentration  at  the  substrate  vs.  pressure  at  fixed  a  and  Too,  assuming 
chemical  equilibrium  outside  the  boundary  layer.  The  dashed  curves  show  the  results  if 
boundary  layer  chemistry  is  neglected. 


FIG.  3.  Predicted  scaling  of  process  efficiency  f]g  with  flow  Mach  number. 
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